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Abstract 



The microscopic mechanisms for universality of Casimir forces between 
macroscopic conductors are displayed in a model of classical charged fluids. 
The model consists of two slabs in empty space at distance d containing clas- 
sical charged particles in thermal equilibrium (plasma, electrolyte). A direct 
computation of the average force per unit surface yields, at large distance, 
the usual form of the Casimir force in the classical limit (up to a factor 2 
due to the fact that the model does not incorporate the magnetic part of the 
force). Universality originates from perfect screening sum rules obeyed by the 
microscopic charge correlations in conductors. If one of the slabs is replaced 
by a macroscopic dielectric medium, the result of Lifshitz theory for the force 
is retrieved. The techniques used are Mayer expansions and integral equations 
for charged fluids. 
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1 Introduction 



It is well known that the fluctuation-induced forces between macroscopic conductors 
have an universal character : they only depend on the shapes of the bodies, but not 
on their material constitution. This observation originates from the celebrated paper 
of H. B. G. Casimir calculating the force between two parallel metallic plates due to 
the fluctuations of the quantum electromagnetic field in vacuum at zero temperature. 
The literature produced since then is so vast that we will only quote in the sequel a 
number of papers relevant to our purpose. Many references can be found for instances 
in the books and reviews [1] [2] [3] [4]. 

Concerning the extension of Casimir's result to non zero temperature T, Balian 
and Duplantier provide the general form of the free energy in presence of ideal 
conductors of arbitrary shapes [5] . The theory of Lifshitz and coworkers generalizes 
the calculations to dielectric bodies caracterized by their dielectric fonctions [6], [7] 
(see also [9]). The ideal conductor situation can be recovered from the latter theory 
by letting the dielectric constants tend to infinity. From these studies one can obtain 
the asymptotic behaviour of the attractive force between two planar conductors at 
distance d at high temperature (or equivalently at large separation d) 1 

k B T((3) 

where ( is the Riemann (^-function. In this regime the force is exclusively due to 
thermal fluctuations and the result may be called classical since it does not depend 
on Planck's constant. In the above mentioned theories the conductors are treated at 
the level of macroscopic physics. In fact they are represented by surfaces, called ideal 
conductors, on which the electromagnetic field has to satisfy the metallic boundary 
conditions. The purpose of this work is to gain an understanding of the microscopic 
mechanisms in the conductor that lead to the universality of the force (1). 

To this end we analyse a simple model where the conductors are described in 
fully microscopic terms. The conductors consist of two slabs at distance d contain- 
ing fluids of classical charges (e.g. classical electrolytes or plasmas). The slabs are 
globally neutral but their material composition (charges and masses of the particles) 
can be different. The space external to the slabs is empty. The system of the two 
slabs, considered as a whole, is at thermal equilibrium with a Gibbs weight that 
includes pairwise interactions between all the particles, consisting of Coulomb po- 
tentials plus short-range repulsions. In this setting we present an exact computation 
of the asymptotic behaviour of the average force per unit surface between the two 
(infinitely-thick) slabs giving 

fc B TC(3) 

</>~-g^3-> (2) 



1 The relevant dimensionless parameter is d ksT/hc, c is the speed of light, H the Planck constant, 
ks the Boltzmann constant. 
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It is also checked that (2) still holds for slabs of finite thickness (appendix A). 

One may notice that the usual approaches consider the fluctuating electromag- 
netic fields as primary objects. The thermal fluctuations of these fields originate 
from the fact that they are in equilibrium with the matter constituting the con- 
ductors, but as a consequence of universality, the microscopic degrees of freedom 
of the charges in the conductors do not need to be explicitly incorporated in the 
description. Here we adopt another point of view : we start from the thermal con- 
figurational fluctuations of the charges to provide a direct calculation of the average 
force without recourse to the field statistics. Then the origin of universality can be 
traced back to the specific sum rules obeyed by the correlations of Coulombic matter 
[10]. 

Universal properties of a variety of classical models of conductor have been stud- 
ied in [11] [12]. In [11] the authors consider a statistical mechanical system of charges 
confined to a plane at distance d of another ideal planar conductor and establish the 
result (2). In [12], they show that replacing the above ideal conductor by fluctuating 
charges does not alter (2). 

The value (1) arising from the electromagnetic field fluctuations calculations is 
twice larger than that obtained in the purely electrostatic models considered here, 
as well as in [11] [12]. This point has been the subject of several discussions in 
the literature, in particular in [8] [9]. In [8] Schwinger performs a calculation of 
the Casimir force mediated by scalar photons (corresponding to the sole electric 
degree of freedom of an electromagnetic wave) leading to the result (2). In [9] the 
authors show that taking the magnetic degree of freedom of the field property into 
account multiplies the expression (2) by 2. In its very formulation our model does 
not include the magnetic part of the Lorentz force induced by fluctuating currents 
in the conductors, whose effect has the same magnitude as that of the Coulomb 
force. Although such purely electrostatic models of conductors do not account for 
the physically correct value of the force at large distance, they already nicely reveal 
the microscopic mechanisms occurring in conductors that guarantee its universality. 

The calculation of the force requires the knowledge of the charge correlation 
function across the two slabs separation, which is the main object of our study. 
At large separation it remarkably factorizes into three parts. There is a first factor 
independent of the slabs' material constitution and two other factors, each solely 
associated to one of the conductors. More precisely, the latter factors involve the 
charge density of the screening cloud induced by a charge located at the boundary 
of a single conductor in empty space. Then the universality of the force results from 
the perfect screening sum rule that holds in any conducting phase. 

In this work we use the technique of Mayer expansion and integral equations for 
charged fluids. In section 2, we specify the system under study and express the force 
by unit surface between two infinitely extended slabs in terms of the microscopic 
charge correlation between them, taking the existence of the thermodynamic limit 
for granted. 

The general formalism used is recalled in section 3 : the charge correlation func- 
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tion is written in terms of the Ursell function subjected to a Mayer expansion. The 
prototype graphs entering in this expansion involve screened Coulomb bonds result- 
ing from chain summations (Debye-Huckel mean-field potential) and density weights 
at vertices [13] [14]. The weights are the exact inhomogeneous densities that have to 
be self-consistently determined from the first BGY equation. We do not treat here 
the full self-consistent problem because it turns out that the detailed structure of 
density profiles is not needed (see [15] [16] for a thorough study of density profiles 
near boundaries). We only have to introduce weak and plausible assumptions on the 
convergence of the profiles to their bulk value. 

It is shown in section 4 that the asymptotic value (2) of the force is already 
obtained at the level of the Debye-Huckel theory. The main tool is the explicit 
form of the mean-field potential for piecewise-flat density profiles, related to the 
potential for structured profiles by an integral equation. The latter equation is shown 
to have a convergent perturbative solution in the weak-coupling regime. At large 
slab separation, the Debye-Huckel potential factorizes into potentials pertaining to 
individual plasmas obeying electroneutrality sum rules. 

We establish in section 5 that the theory beyond mean-field does not provide any 
additional contribution to the asymptotics (2). This is first done non-perturbatively 
with the help of integral equations corresponding to an appropriate dressing transfor- 
mation of the Ursell function and under the mild assumption of integrable clustering. 
Finally the result is recovered once again by selecting and resumming the contribu- 
tion of dominant Mayer graphs to the full charge correlation at large separation. 

We also treat in section 6 a variant situation where one of the slabs is replaced by 
an ideal macroscopic dielectric medium at distance d (namely generating images of 
the plasma's fluctuating charges). Using the Green function of the Poisson equation 
with appropriate dielectric boundary conditions, the Lifshitz result for the mean 
force is retrieved and reduces to (2) as the dielectric constant tends to infinity. It is 
interesting to observe that it is sufficient for the fluctuations to occur only in one of 
the bodies to generate the same asymptotic behaviour. 

We will come back to the inclusion of magnetic forces and to quantum models 
in future works. 

2 Description of the model 

We consider two plasmas A and B of classical point charges confined to two planar 
slabs A A (L,a) and A B (L,b) in three-dimensional space. The slabs have thickness a 
and b, surface L 2 , and are separated by a distance d : 



The plasma A (B) is made of charges e a (ep) of species a G Sa (P G S b ) where 
Sa and S B are index sets for the species in Aa(L, a) and A B (L,b) respectively. We 



A A (L, a) : = {r = (x, y) | x G [-a, 0], y G 
A B (L, b) := {r = (x, y) | x G [d, d+ b], y 



[-if] 2 } 



(3) 
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assume both plasmas to be globally neutral, i.e. carrying no net charge, 



5>«-=Z>A = ° (4) 

a b 

where ^ a extends on all particles in A A (L, a) (A B (L, b)). For a particle located 
at r we will use the generic notation (7, r) where 7 G S A if r € A A (L, a) and 7 G S B 
if r G A B (L, b). The space external to the slabs is supposed to have no electrical 
properties, its dielectric constant being taken equal to that of vacuum. The charges 
are confined in the slabs by hard walls that merely limit the available configuration 
space to the regions (3). 

All particles interact via the two-body potential 

V(l, i, l r - r 'l) = e 7 e y i;(r - r) + w SR (7, V, |r - r'|), (5) 

where v(r — r') = 1/jr — r'| is the Coulomb potential (in Gaussian units) and 
^sr(7, V, r— r') is a short-range repulsive potential to prevent the collapse of opposite 
charges and guarantee the thermodynamic stability of the system. 

The total potential energy U consists in the sum of all pairwise interactions, 
separated into three contributions according to whether they take place between 
two particles of A, of B, or between a particle of A and a particle of B : 

U = U A + U B + U AB . (6) 

On the microscopic level, the force between configurations of charges in the two 
plasmas is the sum of all pairwise forces exerted by the particles of B on the particles 
of A : 



a b 



e aa ep b —^— — - + F SR (a a , A,, r a - r b ) 



r a GA A (£,a), r h eA B (L,b) (7) 



and Fsr is the force associated to the short-range potential i>sr- For simplicity we 
assume that the range of i>sr is finite so that F S R(o;a, A>, r a _ r b) vanishes as soon 
as d is large enough, and we will omit it in the following. 

Both plasmas are supposed to be in thermal equilibrium at the same temper- 
ature T. The statistical average (•••).£, is defined in terms of the Gibbs weight 
exp(— /3U), j3 = {k B T)~ x , associated with the total energy (6). There is no need to 
explicitly specify the ensemble used here (canonical or grand canonical) provided 
that the global neutrality constraint (4) is taken into account. The average parti- 
cle densities p^(7r) are expressed as averages of the microscopic particle densities 
p(7r) = J2i ^77i^( r — r i) where the sum runs over all particles 

p L ( 7 r) = (p( 7 r)) L . (8) 



5 



We keep the index L to remember that averages are taken for the finite- volume slabs 
(3). Hence expressing the sums in (7) as integrals on particle densities p(jt), the 
average force reads 

( F )l= /dr fdr> -Lz£L CL (r,rO (9) 

Ja a (l) Jk B (L) |r - r I 

where Ci(r, r') is the two-point charge correlation function 

c L (r, r') = (c(r)c(r')>L = ^ e 7 e y p L (7 r, 7 r'), c(r) = ^ e 7P(7 r). (10) 
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We now consider the average force by unit surface between two infinitely extended 
slabs at distance d by letting their transverse dimension L tend to infinity. We 
assume that the plasma phases are in fluid states homogeneous and isotropic in the 
y directions, namely the charge correlation has an infinite-volume limit of the form 

lim c L (r,r') = (c(r)c(r')) = c(x, x', |y - y'|). (11) 

L^oo 

For symmetry reasons, {F) L has no transverse component and is directed along the 
x axis perpendicular to the plates. We therefore consider the x-component of the 
force per unit surface 

(/) := lim = lim -J- fdy ( fdx f dx' fdy' ^ — ^ c L (x, y, x', y') 

l^oo L 2 l^oo L 2 J L 2 \J_ a J d J L 2 |r - r'| 3 

-0 rd+b r X — x' 

dx I dx' dy- — —— c (x,x',\y\). (12) 

J [{x - x'y + |y| 2 J 



a J d 



The last line results from the y translational invariance of the integrand in the limit 
L — > 00. We do not justify the existence of the limit here (which depends on an 
uniform control of cl(x, y, x' , y') as |y — y'| — > 00), but it will be clear from the 
subsequent calculations that (12) is a well defined quantity, at least in the weak- 
coupling regime. 

Formula (12) remains valid if one replaces c(x,x', |y|) by the truncated charge- 
charge correlation function 

S(x,x',y) = (c(r)c(r'))-(c(r))(c(r')), r = (x,y), r' = (x',0) (13) 

with c(r) the microscopic charge density as in (10). Indeed, the y- Fourier transform 
of the Coulomb force reads 



■ [ y ^' y [{x -x>) 2 l\y\^ = 2n sign(x - x ' )e ~ klx ~ x ' 1 (14) 
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and reduces to — 2ir when k = and x < x'. This implies that the charge density 
profile (c(r)) = c(x) does not contribute to the force because of the global neutrality 
of both plasmas 

/•0 pd+b 
dxc(x) = / dxc(x) = 0. (15) 
-a Jd 

To take full advantage of the translational invariance in the y direction we rep- 
resent the y-integral in (12) in Fourier space : 

(/) = -i- Jdx J^dx' Jdke- k \ x - x '\S(x,x',k) (16) 

where k = |k| and S(x,x',k) = fdke~ tk ' y S(x,x',y). The dependence of (/) = 
(/) (d) on the separation d between the two slabs occurs in the integration limits in 
(16) as well as in the charge correlation function S(x,x',k). The d dependence of 
the correlations between the two slabs A and B originates itself from the Coulomb 
interaction term Uab occurring in the total Gibbs thermal weight. The object of the 
next sections is to determine the asymptotic behaviour of (/) (d) as d — > oo. 



3 Mayer series for inhomogeneous charged fluids 

We briefly summarise the methods that we use to calculate the charge-charge cor- 
relation function of our system. Let us consider a general charged fluid in presence 
of spatial inhomogeneities caused by an external potential f cxt (7r), e.g. wall po- 
tentials confining the system in some region of space. Hard walls without electrical 
properties (infinite potentials) can be implemented by simply declaring that the 
density vanishes in the forbidden regions. 

It is well-known (e.g. [17]) that the two point Ursell function, related to the 
densities p(i),p(j) and the two-particle correlation p(i,j) 

h(iJ):=-$tfL-l, (17) 

can be expanded in a formal power series of the densities by means of Mayer graphs. 
The basic Mayer bonds are 

f(i,j) = e-^-l (18) 

where V(i,j) is the potential (5) and the weights at vertices are the densities p(i). 
Here i is a shorthand notation for the point (7* r«) in configuration space, and in- 
tegration on configurations ^ f d?i includes the summation on particle species. 
Diagrams have two root points % and j and m internal points which have to be in- 
tegrated over. Each pair of points is linked by at most one /-bond and there are no 
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articulation points 2 . Because of the long-range of Coulomb interaction, the integrals 
occurring in every diagram diverge in the thermodynamic limit. It is therefore neces- 
sary to introduce the screened mean- field potential r 3 -) as usual by resumming 
the chains built with pure Coulombic interaction bonds — f3e li e lj v (r^ — r^). Then 
replacing the bare Coulomb potential by the screened potential leads to a reorgani- 
sation of the diagrammatic expansion of the Ursell function resulting in the formula 
[14] 

h(7r,7'r') = ^— ^ /dri • • • dr m p( 7l ri) • • • p( 7m r m ) JJ 
n n 7 i,..., 7m ^ {i,j}en 

(19) 

The first sum runs over all unlabelled topologically different connected diagrams II 
(called prototype graphs) with two root points (71") and (7'r') and m integrated 
internal points with density weights (m ranges from to 00); Sfi denotes the sym- 
metry number of a diagram II. Each pair of points is linked by at most one bond 
T G {F,F R } and there are no articulation points. Moreover, convolutions of F 
bonds are forbidden to avoid multiple counting of original Mayer graphs. The two 
possible bonds read 

F{i,j) = -/3e Ji e y ^(T i ,T j ) (20) 

F R (i,j) = exp[-/3e 7i e 7 .$(r;, r,) - pvsn^i, 7j, |r s - r^-j)] - 1 + /3e 7i e 7j $(r;, Tj). 

(21) 

These bonds are obtained in terms of the Debye-Hiickel screened potential $, which 
is symmetric and defined as the solution of the integral equation 

$(r, r') = v(r - r) - [d^ /t 2 ( ri ) v(r - r : ) $(n, r') = $(r', r) (22) 

4tt J 

or equivalently of the differential equation 

A$(r, r) - /t 2 (r)$(r, r') = -An5(r - r) (23) 
supplemented by suitable boundary conditions. In (22) and (23) 

«(r):= ^47r/3^e 2 p( 7 r)j (24) 

can be interpreted as the local inverse Debye screening length in the inhomogeneous 
system. The bond F R (i,j) includes the short-range contribution and the nonlinear 
Coulombic part of the original Mayer bond. 

2 An articulation point, when removed, splits the diagram into two pieces, at least one of which 
is disconnected from the root points 
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The densities p( 7 r) entering in (19) and (24) have to be determined self-consistently 
from the first equation of the Born- Green- Yvon hierarchy which links the one-point 
and the two-point functions. For charged systems, it takes the form [10] 



Vp(7 r) = - /?e 7 p(7 r) 



W( 7 r) + J At' f^eyp( 7 V)/i(7r,7V) j Vu(r-r') 



-Pj2j dr ' r ' ^ r ') V^sr(7, 7', |r - r'|) 



(25) 



where 



( 7 r ) = ^ cxt ( 7 r) + Jdr' c(r>(r - r') (26) 

is the sum of the external potential and the electrostatic potential caused by the 
inhomogeneous mean charge density c(r') in the system. Hence the Ursell function 
(considered as a functional of the densities through its Mayer expansion (19)) to- 
gether with (25) form a closed set of equations whose solution determines in principle 
the exact densities and two-particle correlations. The differential equation (25) has 
still to be supplemented with appropriate boundary conditions. For instance if the 
system is asymptotically uniform in some directions, one can fix the corresponding 
asymptotic bulk densities. 

Finally, the charge-charge correlation function (13) is related to the Ursell func- 
tion by 

S(r, r') = <WA>(7 r)p( 7 ' r')M7 r, i r') + 8{v - r') £ e>( 7 r) (27) 

7:7' 7 

The second term in the r.h.s. of (27) is the contribution of coincident points. 

We shall use the above formalism to calculate (27) as a function of the distance 
d for two infinitely-thick plasmas (i.e a — > oo and b — > oo in (16)) and with hard 
walls at x = and x — d. In this situation we take \I> cxt ( 7 r) = for x < and 
x > d, and impose 

p( r/ x) = 0, 0<x<d, 
lim p( 7 x)=p j47 , lim p( 7 x) = p Bl (28) 

where p^ 7 and p^ 7 are the bulk particle densities of plasmas A and B. The con- 
tribution of coincident points does not enter into the force (16) since r and r' are 
always at least separated by the distance d. Therefore, (16) reads 

/•O roc poo 
da; / da;' / dA; ke'^-^ ^ e 7 e y p( 7 x)p{i x')h(>y x, i x', k) (29) 
-oo Jd JO , 

with h^x, 7' a;', k) the y- Fourier transform of the Ursell function. 
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4 Debye-Hiickel theory 

In this section we show that the simplest contribution to r, 7' r') given by the 
sole bond F, namely, 

/i DH ( 7 r, 7' r') = -/3e 7 e y $(r, r) (30) 

already leads to the asymptotic value (2) of the force. For this we have to find 
the screened potential by solving (23) (written in Fourier form) with the boundary 
conditions imposed by the slab geometries 



<9 2 

-k 2 -K 2 (x) 



dx 2 



<S>(x,x',k) = -4n5(x - x'), k(x) = 0, < x < d, (31) 



with $(2;, x',k) = Jdy e _lk ' y <E»(a;, x', y). The boundary conditions are <&(x, x', k) 
and d$(x,x' ,\c)/dx continuous at x = and x = d, and lim. I ._ >±00 $(x, a;', k) = 0. 
The density profiles entering in k 2 (x) by (24) are not known (since they have to 
be determined by self-consistency from (25)), but we will not need their explicit 
form in the sequel 3 . We only need to assume that their difference to bulk value is 
integrable : 



p{ix) - pB~ t = O , x -> 00, e>0. (32) 

Integrating (31) on x leads to 

da; — —Q(x, x' , k) = 1 - — /dx a/, k). (33) 

-OO ^ «/ 

In particular, for k = 



da;— y -±${x,x',k = 0) = 1 (34) 

47T 



which is nothing else than the electroneutrality sum rule for the charge-charge cor- 
relation (13), (27) within the Debye regime (30) [10]: 

[dx [dyS DH (x,x',y) = (35) 



3 A mean- field approximation to the densities could be obtained by replacing h in (25) by ft, DH . 
We are not doing so here but deal throughout with the exact densities. 
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To solve (31) we first consider the simpler problem with piecewise-flat densities 
PA-y and pb<j in each plasma 



<9 2 

k 2 -R 2 {x) 



... , ip(x, x' , k) = — 4tt5(x — x) (36) 

ox 

k(x) = ka, x < 0, R(x) = 0, < x < d, R(x) = Kb, x > d 

where 

K A = I 4ltf3 ^ 6 l PAa ) , K B = I 4:71(3 ^2 e l Pbp ) ( 37 ) 
V aeS A J \ I3&S B ) 

are the bulk inverse screening lengths. The boundary conditions are the same as for 
(31). Denoting C the linear operator acting on $ on the l.h.s. of (31) and by L the one 
acting similarly on ip, one has £$(x) — Cip(x) = £($ — <p)(x) — u(x)$(x) = 0, where 
u(x) = k, 2 (x) — R 2 {x) represents the deviation of the density profiles to their bulk 
limiting values. Since —(p/4n is the Green function of C, it follows that $(x,x',k) 
and ip(x, x', k) are related by the integral equation 

$(x,x',k) = ip(x,x',k) - -j- fdsu(s) tp(x, s,k) $(s,x',k) (38) 

4lX J 

which expresses $(x,x',k) as a perturbation of <p(x,x',k) by the inhomogeneity 
u(x) of the plasmas' density profiles. 

Solving (36) piecewise and connecting the solutions together yields 4 

{ip AA {x,x,k), x,x'<0, 
(Pab(x, x' — d, k), x < < d < x', (39) 

<Pbb(x — d, x' — d, k), d<x,x', 



with 

k A \x— x'\ „— kj\\x+x'\ 
V AA(x,x',k)=2n- + 27T- fr.-^^+^-q-,,-/.-)^-/..),. 



e , o e (k A ~k)(k B +k)c kd -(k A +k)(k B -k)c- kd 

fc A ^ fc A (k A +k){k B +k)e kd -{k A -k){k B -k)e- kd > 



(40) 

ip AB (x, x', k) = jj- - - ^ kd - jj- _ -jy- _ ^_ kd , (41) 

k A = \jk 2 + K 2 A , k B = y/k 2 + 4- (42) 



4 The functions <p ab, Vbb (depending on d) refer to the system of the two plasmas under 
mutual influence with the allocation of particles in plasma B measured by their distance from the 
boundary at d (i.e. from to +oo). In the sequel, the quantities ip A , <p% (independent of d) refer 
similarly to the single semi- infinite plasma A and B. 
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The function ipbb(x, x', k) is obtained by interchanging the indices A and B in 
(40). Notice that ip(x, x', k) = </?(V, k) and is invariant under the symmetry x <-> 
d — x, x' <-> d — rr', A <-> £>. 

We discuss a few properties of this solution. The first term in the r.h.s of (40) 
corresponds to the bulk Debye-Hiickel potential whereas the second term is the 
modification due to the finite boundaries of both plasmas A and B. As d — > oo, 
(Paa(x,x',\s) reduces to the well-known Debye-Hiickel potential (p A ( x , x', k) of a 
single semi-infinite plasma in the region x < (see formula (24) in [18] and [19]) 

e ~k A \x-x'\ _ ^ e -k A \x+x'\ 

lim <£aa{x, x', k) = 2vr h 2vr- — = <p° A (x, x', k) (43) 

uniformly with respect to k. One observes that ip(x,x',\s) is an even, infinitely 
different iable function of |k|, implying that tp(x,x',y) decays along walls directions 
faster than any inverse power of y. This is to be contrasted with the small k behaviour 
of the function (43) which has a non analytic |k| term leading to the algebraic decay 
y~ 3 along the wall [19]; see also [10], sect. III.C.2. 

The function <pab{%, x ', k) (41) describes the correlation between the two plas- 
mas. In terms of the scaled dimensionless variable q = k d, it has the simple factor- 
ized asymptotic behaviour 

<Pab(x,x',*) ~ ^-e-^le^l 

= ^ 47r J nh ^(x, 0,0)^(0, ^,0), d^oo (44) 

Finally ip(x, x', k) obeys the following bound uniformly with respect to k and d 
(appendix B) 

4:71 

< ip(x,x',k) < ip > (x,x') < — , nd>\ (45) 

K 

The function ^(x, x') is defined as in (39) with <paa, ¥ab and ip BB replaced by 

<P$a(x, =4= &b(x, =4= f (e-^^' 1 + e"^' 1 ) (46) 

<p^u(x, x') = <p> A (x, x') = -^L e —\*\ e —W\^ K ■- min{K A , n B } (47) 

n z d 

The Debye-Hiickel potential x', k) can be obtained by iterating the integral 
equation (38). Convergence can be established in the weak-coupling regime : 

Lemma (see proof in appendix B) 

Let 



r :- 



^Jdx\u(x)\ = ^Jdx\K 2 (x) -R 2 (x)\ (48) 
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Then for r < 1 (38) has a solution with the bound 

\^ XiX \k)\<-^^{x,x'). (49) 

As in (39) we distinguish various contributions according to the location of the 
arguments x, x' of $(x, x',k) by setting (see footnote 4) 

{$ AA (x,x',k) , x,x' < 0, 

$ AB (x,x' - d, k) , x < < d < x', (50) 

^_bs(^ — d, x' — d, k) , d<x,x', 

( \ \ Pa(ix) x < 0, 

P(7^) = < , ,v , (51) 

yp B (^,x -d) , x > d. 

The quantities k AjB (x), ua,b(x) are defined in the same way. Then the integral 
equation (38) splits into 

1 f° 

$aa(x, x', k) = ip AA (x,x',k) - — <Isua(s)<paa(x, s, k)$ AA (s, x', k) 

4tt J_ 00 

*<X> 

dsM B (s)v? AB (x, s, k)$ BA (s, x', k). (52) 



1 

47T 



i r 

$As(^,^',k) = ^As(^,^',k) - — dsU A (s)lf AA (x, S, k)$ AB (s,x',k) 

1 Z" 00 

-— / dsu B (s)(p AB (x,s, k)$ BB (s,x',k). (53) 
47r Jo 

The density profiles depend on d because of the mutual Coulomb interactions 
between the two plasmas. We shall examine the asymptotic behaviour of $(x, x' ', k) 
as d — > oo under the assumption that these density profiles are uniformly bounded 
with respect to d and tend to those of single semi-infinite plasmas i.e. 

lim p A {^x) = p° A (ix), x < 0, lim p B {^ x) = p° B (-fx), x > (54) 

a— >oo d^oo 

We denote by k ab (x), u ab (x) the analogous quantities for the single semi-infinite 
plasmas. Then one concludes from (52) that uniformly in k 

lim &aa(x : x , k) = $° {x,x',k), x,x'<0 (55) 

d— >oo 

where & A {x, x', k) is the Debye-Hiickel potential of a semi-infinite plasma in the re- 
gion x < determined in terms of the corresponding flat profile potential (p A (x, x', k) 
(43) by 

i r° 

<5> A (x,x',k)=p A (x,x',k)- — dsu A (s)(p A (x,s,k)<$> A (s,x',k), x,x'<0 

(56) 
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Indeed, in view of the limits (43), (54) and using dominated convergence with the 
bounds (45), (49) the integral equation (52) reduces to (56) in the limit d — > oo. 
One has likewise 



lim $ BB (x,x',k) = $° B (x,x',k), x, x > 



d^oo 



(57) 



where $ B (x, x',k) is the Debye-Hiickel potential of a semi-infinite plasma in the 
region x > 0. 

We come now to the correlation &ab(x, x', f ) which is expected to decay as d^ 1 
at large separation of the two plasmas. To see this it is useful to write (53) in an 
alternative form such that ifiAB appears explicitly in each term of the equation : 

$ab(x, x', k) = (Pab(x, x', k) 

1 f° 1 f°° 

— ds u A (s) ^AA(x,s,k)(f A B(s, x' , k) - — / ds u B (s) <pab(x, s, k) f BB (s, X 



4ir 



47T 



+ ih) I ds i / ds 2 u A (si)u B (s 2 ) ®aa(x, si, k) v?ab(si, s 2 , k) $bb(s 2 , x', k). 

-oo JO 



(58) 



Here <&aa(x, x', k) verifies the equation 



$AA{x,x',k) = Lp A A(x,x',k) - — dsu A (s)(p A A(x,s,k) $aa(s,x', k), (59) 

471" J-oo 

which is the equation (52) with the last term omitted. Equation (58) is obtained by 
iterating (53) and resumming the <£aa chains to §aa, or by verifying that it satisfies 
the basic differential equation (31). By the same arguments that led to (55), it is 
clear that also tends to the potential $^ of the semi-infinite plasma. 

Introducing the limits (44), (54), (55), (57) in (58) and using again dominated 
convergence provided by the bounds (47), (46), (49) we find that 



lim d &ab(x, x\ 

d—*oo 



d' 



Q 



Air sinh q 



x 



x 



V ° A (x, 0, 0) - — jds u A (s)K(x, s, 0)<p° A (s, 0, 0) 



oo 
oo 



V^(0, x', 0) - — / dsu° B (s)<p° B (0, s, 0)$° B ( S , x', 0) 



q 



47r sinh q 



4tt Jo 

$°(x,0, 0)^(0, x',0) 



(60) 



As in (44), the limit factorizes into the product of Debye-Hiickel potentials for single 
semi-infinite plasmas evaluated with one point on the boundary. The last line of (60) 
follows from (56), the corresponding equation for $ B (x, x', 0), and the fact that these 
functions are symmetric in x, x' . 
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With this result we can determine the leading term in the asymptotic behaviour 
of the force (16) in the Debye-Huckel regime. From (29) and (30), one has 

(/>DH {d) = 1 ft l?[ Akk ^"" ^^*«"' k) 

To obtain the second line we have set k d — q, shifted the x'-integration by —d, and 
introduced the notation (50), (51). As d — > oo, one can use (60) and the bounds (46), 
(49) to conclude again by dominated convergence that 

^oo w ' y ' 871,3 J e- - e-> Vi_o» 4tt m ') 

"OO 



Indeed, because of the charge sum rules (34) for the semi-infinite plasmas, the two 
parenthesis are equal to 1, whereas the q integral yields the value £(3) where ( is 
the Riemann (^-function. 



5 Contributions of the other graphs 

In this section we show that the single Debye-Huckel bond F saturates the asymp- 
totic behaviour of the force i.e. taking into account the full set of other diagrams 
does not modify the result (62). For this we use the method of "dressing of the root 
points" that has been introduced in [20] to analyse the decay of the quantum trun- 
cated charge correlation function. Having singled out the contribution of the single 
resummed bond F, all remaining diagrams constituting h{^f r, 7' r') can be classified 
into four types, depending on whether their root points are linked to the rest of the 
diagram by a single bond F or not 5 . We thus write their sum in the form 

h R := h - F = h cc + h cn + h nc + h Dn , (63) 

where h cc stands for the contribution of all graphs that do begin and do end with 
an F bond (with anything in between), h cn for the contribution of those that do 
begin but do not end with an F link, and so on. The latter quantities are obviously 
related to h im by the following integral equations (notations are as in sect. 3) 

h cn (a,b) := JdlF(a,l)p(l)h nn (l,b) 

h cc (a,b) := Jdl Jd2F(a,l)p(l)h nn (l,2)p(2)F(2,b) (64) 



5 A point in a prototype diagram which is linked to the rest of the diagram by exactly one F 
bond is called a Coulomb point in [20]. 
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and analogously for h nc . Using these representations in (63) together with (20) and 
building the charge-charge correlation corresponding to h R according to (27) yields 

S R (x,x',k) = J2 [ dx iYl [ dx2 (s(x-xi)-^^Hx, Xl ,k)\ 



71 72 



x e 7iP(7i xi) ^ nn (7i x i,l2 x 2 , k) e 72 p(7 2 x 2 ) ( 5(x 2 - x') -^—§(x 2 , x , k) j . 

(65) 

The function h nn {ji x±, 72 x 2 , k) embodies a resummed contribution, not explicitly 
known at this point, of higher order graphs. The only assumption needed on this 
function in the sequel is integrable clustering uniformly with respect to d 

"OO 

dxi /i nn (7i xi, 72 x 2 , k) < oo (66) 

poo 

dxi / dx 2 /i nn (7i ^1,72X2, k) < cx) (67) 

-00 J d 

As a consequence of the bounds (45), (49) the condition (67) obviously holds for the 
Debye potential $, and it is expected to hold for the Ursell function on the ground 
that as x\ — > —00 (x 2 — > 00) fo nn (7i xi, 72 £2, k) has a fast decay in the bulk part of 
plasma A (plasma B). Note that integrating (65) on x (or x') at k = proves the 
validity of the charge sum rule for the exact charge-correlation function (see (35)) 

/'OO f 
dx dy S(x, x', y) = 0, S(x, x', y) = S DH (x, x', y) + S R (x, x\ y) (68) 
-00 J 

Proceeding as in (61) the contribution of S R (x, x', y) to the average force can be 
written in the form (permuting the x, x' and Xi,x 2 integrals) 

(f} R (d) = -^JdqqJdx 1 Jdx 2 H l {x 1 ^) 
x e J1 e l2 p{-f 1 x 1 )p(-f 2 x 2 )h nn (-f 1 x l , -y 2 x 2 ,^,d)H 2 (x 2 ,^) (69) 

71:72 

where 

H 1 (x 1 , 2) = jdx (s(x - x ± ) - ^-®(x, x u §)) 

H 2 (x 2 , 2) = jdx' (6(x 2 - x>) - ^1^( X2 , X ', *)j e"^' (70) 

The behaviour of (/) R (d) as d — > 00 is determined by that of the functions Hi and 
H 2 , because /i nn (7i xi, 72 x 2 , ^, d) does not vanish in the limit when the variables 
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x\ and X2 are both located in the same plasma, but tends to the corresponding 
functions associated with a single semi-infinite plasma. Both Hi and H 2 are 0(1/ d) 
so that {f) R (d) = C(l/o? 4 ) does not contribute to the asymptotic behaviour of the 
force. More precisely, integrating (31) on x gives for x\ < 

y° N k 2 (x)^ A « 2 (x)_ A; 2 /" _ 

dx [8{x - xi) -^$(x, xi,k) )= / dx <E>(x, x i; k) + — dx$(x,: 

J-oo V 47F J Jd 471 47T J 

dx x 1 ,k) + —jdx $(x, x u k) (71) 



= dx 
Jo 

implying with (46), (49) 

r 

d7 1 ~ Vrf 2 



H 1 (,- 1 .H) = 0(1) • ' ( "F ) - <0 (721 



For xi > d one has 

e -« B (xi-d) ( 73 ) 



^(xx, 1) = - fdxe«*' d ^*(x, Xl , a) = O Q) 



In the same way 

ff 2 (x 2 ,s}) = 0^) e -^l, x 2 <0, 

^ 2 (^l) = o(^)+o(^), x 2 >ci (74) 

(the factor e _<? comes from e~ qx / d < e~ q for a; > d in (70)). Inserting these estimates 
in the four integration domains determined in (69) by xi, x 2 < 0, x±, x 2 > d 
together with the integrability assumptions (66), (67) on h nn leads to the result 

Jim d 3 (f) (d) = Jim d* (/) DH (d) + Jim d 3 (f) R (d) = ^ (75) 

d^oo d^oo d^oo oTtfJ 

To conclude this section we present an alternative derivation of the result (75) 
by selecting the class of diagrams that give the dominant contribution to the Ursell 
function as d — > oo. For this we decompose the bond F(p(x, "f'x 1 , k) (in Fourier 
representation) into the sum of four terms according to the location of the arguments 
x, x 1 



F = Faa + Fab + F B a + F B b 

F(~fx,Y x',k) , x,x' < 
, otherwise 



F AA (-fx,j'x',k) 



- , , . F(-fx,-f',x' + d,k) , x<0,x'>0 
F AB {^x,-f x ,k) = < (76) 

, otherwise 
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with Fba and F BB defined likewise and the similar decomposition for F R (re-integrals 
in plasma B from now on run in the interval [0, oo), see footnote 4). The set of 
prototype graphs is then expanded in a larger set of graphs defined in terms of these 
bonds. It follows from (55) that Faa and F AA bonds have limits F A and F A R as 
d — > oo where F A and F A R are the bonds corresponding to the semi- infinite plasma 
alone and likewise for BB bonds. 

It is shown in appendix C that the dominant part of the Ursell function h AB {p( a T a, lb 
as d — > oo is constituted be the set of graphs that have exactly one F A b bond. This 
class is obtained by linking the extremity 71 ri of Fab (71 72 r 2 ) to the root point 
7 a r a of h AB [^ a i"a, lb *b) i n plasma Aa by all possible subgraphs comprising only AA 
bonds (otherwise one would introduce additional AB bonds), taking into account 
the excluded convolution rule for F bonds. In the same way the other extremity 
72 r 2 of Fab (71 r i, I2 ^2) is linked to the root point % r b in plasma A B by all possible 
subgraphs made of BB bonds. One finds in this way 

h AB (a, b, k) ~ Jdl Jd2 [5(a, 1) + (h n A n A (a, 1, k) + h^a, 1, k)) p A (l)] 

x F AB (1, 2, k) [5(2, b) + (h n B n B (2, b, k) + hf B (2, b, k)) p B {2)] (77) 

Here a = (i a x a ), 1 = (aiXi), 2 = (P2X2), b = (75, Xb), and the integration J dl — 

J-00 mns 011 Plasma A and J d2 — Ylifo Jo°° ^ x 2 on plasma F. We have also 
used that the convolution of translation invariant functions in the y-direction is the 
product of their Fourier transforms. The functions h c AA and h n AA are defined as in 
(63) in terms of AA bonds (similarly for h c BB and h B n B in terms of BB bonds). One 
can write (60) in the form 

FA n 2;3 x g F°(l,a ,0)F°(& ,2,0) ^ ^ ^ 

' ' d A^f3d sinh g e ao ep 

where ao = (ao 0) indexes a charge e ao located at the boundary x ao = of and 
bo indexes a charge at the boundary of A B . Taking also into account that the 
functions h c AA and h n AA approach the corresponding values (h° A ) cn and (h A ) nn of a 
single semi-infinite plasma, one finds that the leading term ~ 1/d of Kab (l x ,l' x ' 
factorizes as 

h (a b 3 ) ~ - G A (a,ao) G° B (b ,b) ^ 

' d 47T/9(isinhg e ao ep 



with 



G A (a,a )=F° A (a,a ) + jdl [(h A )™(a, 1) + (^) cn (a, 1)] p A (l)F° A (l, a ) 

= (F° + + ) CC ) (a, ao) (80) 



In all functions are evaluated for the single semi-infinite plasma A at k = and 
k has been omitted from the notation. The expression for G° B is built in the same 
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way. By the same calculation that led to (62) one finds now from (29) that 

hm d 3 ( f) (d) = ^ ( Jdae aa p A (a)G° A (a,ao) \ ( Jdbe Pb p%(b) G%(b ,b) 

(81) 

It remains to see that both parenthesis are equal to —1 because of the electroneu- 
trality sum rule in semi-infinite plasmas. Indeed, using (63) and (64), one recognizes 
from (80) that 

G° A (a, a ) = h° A (a, a ) - Jdl [F° A (a, l)p A (1) + 5(a, 1)] (h A ) nn (l, a ). (82) 
The contribution to the force of the second term of (82) involves 
Jdae aa p° a (a) [F°(a, l)p° A (1) + o(a, 1)] 

= e Ql/0 (ai xi) (l - ydx a «) 2 (:r a )$° (x , xi, k = 0)^ = (83) 



which vanishes because of the sum rule (34) in the case of a semi-infinite plasma. 
The contribution of the first term of (82) is 

/ da e aa p° A {a) h° A (a, a ) = / dx a e aa p A (^ a x a )h\(a a x a , a 0, k = 0) = -e Qo 

(84) 

The left-hand side is the total charge of the screening cloud induced in the semi- 
infinite plasma A A by the boundary charge e ao , which equals — e ao because of perfect 
screening [10]. By the same considerations the second parenthesis in (81) also equals 
— 1, hence the final result (75). 



6 Plasma in front of a macroscopic dielectric medium 

In this section we investigate the situation where plasma B is replaced by a semi- 
infinite macroscopic medium of homogeneous dielectric constant e. The electrostatic 
potential V(r, r') at r created by a unit charge at r' G is the Green function 
of the Poisson equation with the conditions that the normal component of D(x) = 
e(x)E(a;), (e(x) = e, x > d, e(x) = 1, x < d) and the longitudinal component of 
E(x) are continuous at the interface [21] 

l/|r-r'|+A/|r-r'*|, x < d, A = (1 - e)/(l + e) 
A/|r-r'|, x>d, A = 2/(l + e) 1 ' 
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where r* = (2c? — x, y) is the point symmetric to r with respect to the dielectric 
surface. The case of an ideal grounded conducting plate (e = oo) is formally recovered 
when A = — 1, A = 0. In linear electrostatics, the total energy associated to a 
distribution of charges c(r) external to the dielectric is 

i y"drE(r) • D(r) = X - J^Jdr' c(r)V(r, r')c(r'). (86) 

For a configuration of charges {e Qi , rj} in A^, c(r) = £V e Q , i o T (r — r^), one finds with 
(85) that the total energy can be written as 

U = U A + U AB , (87) 

a N N 



f^A = i r + v S R(a>i, aj, ^ - rj ), /7ab , , , 

where we have omitted the (infinite) self-energies of the particles and added a short- 
range repulsive potential for thermodynamic stability. As in (6), Uab refers to the 
additional energy due to the presence of the dielectric at distance d. The total force 
exerted by the dielectric on the particles of A is obtained by differentiating U with 
respect to d 

F B ^ A = ^U AB = AJ2J2 e ai e aj (88) 

i=i j=i 1 i< 

It corresponds to the sum of all pairwise forces between charges in the plasma A 
and their image-charges Ae Q .. Proceeding as in the derivation leading from (7) to 
(16), the average force along the x direction per unit surface is given by 

=^ftkS^ m ~ i[ ' w ^ s{x ' x '-i ] - (89) 

To obtain the second line we have introduced the Fourier transform (14), used 
x'* = 2d — x', and set k = q/rf. S(x,x',y) is the truncated charge-charge corre- 
lation function of the plasma A defined in terms of the statistical weight exp(— (3U) 
associated to the energy (87). The asymptotic analysis of (89) differs from that of 
the previous section on two points : here the function S provides a contribution 
from coincident points to the integrals and we expect S(x,x', ^,d) to tend towards 
a non-zero limit when d — * oo, namely 

lim S(x,x'S) = S°(x,x',0) (90) 

d—foo 
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with S°(x, x', 0) the charge correlation of the semi-infinite plasma in absence of the 
dielectric. Hence, the leading behaviour of the force as d — > oo 

(/> ~ (/> m ° n - + </> dip - (91) 

comes from the first two terms resulting from the expansion of exp{— q\x + x'\/d} ~ 
1 — q\x + x'\/d in (89). At leading order in (/) dip ' one can replace S(x,x', 3) by 
S°(x,x',0) so that 

roc rO pO _A 

lim d 3 (/) dip ' = -A dq q 2 e~ 2q dx dx' (x + x')S°(x, x', 0) = -. (92) 

d ^°° Jo J-oo J-oo lu7rp 



The first term of the x, x' integrals vanishes because of perfect screening whereas 
the second one equals 1/(4tt/3) as a consequence of dipole sum rule in a semi-infinite 
plasma ([10], formula (3.9), sect. C) 6 . Since lim^oo d 2 (f) mon ' = because of perfect 
screening, one can replace S(x,x', |) by S(x,x',^) — S°(x,x',0) in the monopole 
contribution (/) mon ". It is convenient to further add and subtract S°(x,x', 3) and to 
note that : 

limrf 3 /^ [dx [dx 1 [dqqe- 2q [S°(x,x'Z) - S° {x , x' , 0)}\ (93) 

j}^ ai„ {St fy s °^) = ma- < 94 > 



-A 



This follows from the small k = |k| expansion of the x, x' integral that has a linear 
term ^ ([10], formula (3.24), sect. C). Collecting (92) and (94) in (91) we see that 
the large-d behaviour of (/) is determined by 

-A A f° f° f°° 
(/> = _-- Jdx Jd^ jdqqe- 2 * [S(X,X>, f) - S°(X,X>, *)] + o(rf- 3 ). 

(95) 

One can now proceed as in section 3 with Mayer bonds defined in terms of the Green 
function (85). The Debye-Hiickel equation (31) is supplemented with the boundary 
condition lim^d d$(x, x' , k)/dx = e lim x ^d+ d&(x, x', k)/dx, and k(x) is as before 
for x < 0, and k(x) = 0, x > 0. The solution for piecewise-flat densities is ([15], 
formulae (3.2)-(3.5)) 7 

V {x,x,k)-Zn ^ +271 ^ {kA + k)ekd + A{kA _ k)e . kd - W 



6 Here the sign is opposite to that in [10] because the plasma is located in the x < half space. 

7 Here the minimal distance d between a charge and the dielectric wall plays the role of the 
hard-core diameter in [15]. Notice that we defined A with the opposite sign and that our plasma 
fills the region x < 0. 
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It is convenient to single out the potential <p°(x, x',k) (43) for the semi-infinite 
plasma in the absence of the dielectric and to split 

(p(x, x', k) = (p°(x, x', k) + (Pab(x, x', k), (97) 

where, from (43) and (96), 

Snke~ kA ^e~ kA \ x '\Ae~ kd 
<p AB (x, x', k) = ^- - fc) ^j- - k ^ kd - jj- - k)Ae _ kd] (98) 

One observes that <£ab(x, x', k) has the factorization property analogous to (44) 

lim d(p AB (x, x', 3) = z(q)<p°(x, 0, 0)^(0, x', 0), (99) 

d— >oo 

Z (g) = JL - 100 

The potential <J> corresponding to the exact non-uniform profile is related to tp by 
the integral equation (38). With a reasoning similar to that leading to (58) and (59), 
it can also be split in two parts 

$(x, x', k) = $(x, x', k) + $ab(x, x', k) (101) 

Here <£> verifies eq. (59) with (p° in place of <£aa and tends to the potential $° of the 
semi- infinite plasma in vacuum. <&ab solves eq. (58) with u in place of ua and ub, 
$ in place of &aa, and $ in place of Then using (99) and the fact that both 
$ and <£> tend to $° leads to the asymptotic factorisation of $ab 

lim d$AB(x, x', 3) = z(g)$°(a;» 0, 0)$°(0, x', 0). (102) 

d-^oo 

We are now ready to evaluate the force (95) as d — > oo in the Debye-Hiickel approx- 
imation. From (27), (24) and (30) one has 



s ™ (x , x i, f ) = - l -^M^l Hx y, 3) + % _ ^ (103) 



and the analogous relation for S' 0DH (a:, 3); some care has to be exercised here since 
coincident points do contribute when both x, x' are in the same integration range. 
We subtract and add $(x, x', |) to $(x, x', |) in (103). This gives two contributions 
to the force (95). The first one is 

-A A /\° f° f°? 2o k 2 (x) k 2 (x') r. , q . ~ , q .i 



8ir[3d 3 (3d 2 



-oo JO 

" A + w&&J? q, '* !f £r £ £ 1 * M * ) (104) 



8ir(3d 3 (3d 
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The second one involves the quantity 
1 



47T/9 



[k\x) - ( K °) 2 (x)] 



+ -J ldx> 



i r° , r^«V)~ _ (^) 2 (^) (^Vy , ^ q . 



oc 



|dx' [$(x,x',3)-$ (x,x',3)] =o(Jp) (105) 

This equality follows from the relation (33) for $ (relative to k) and for $° (relative 
to k°) since both potentials satisfy the basic differential equation (31). Thus the 
contribution of (105) to the force is O (Jr). With the factorisation (102) and using 
the sum rule (34) for $°, one finds from (104) and (100) the final result 

J5/ <» = 8^5 + J /d„e-^ W = ^ E { ^J- (106) 

It can be shown along the lines presented in section 5 that the non mean-field part 
of the charge correlation function does not contribute to this asymptotic result. 
According to (101) one splits the bonds F = F + F A b and F R = F R + F% B with 
F%b = exp(F AB - (3v SR ) - 1 - F AB . The bonds F and F R tend to the bonds F° and 
F 0R pertaining to the semi-infinite plasma without dielectric, whereas F A b vanishes 
in the limit. At large separation, F R B ~ (F AB ) 2 vanishes more rapidly than Fab ( see 
appendix A). As in the analysis leading to (79), the leading behaviour of the Ursell 
function comes from graphs having bonds F, F R , a single F A b one, and it takes the 
factorized form (79). The only difference is that both functions G° refer to the same 
plasma A. Then one establishes the validity of (106) as in (81)-(84). 

This result coincides with that of Lifshitz. Indeed, a straightforward generali- 
sation of its asymptotic force as ^> 1 (formula (5.5) in [6]) to the case of two 
different homogeneous dielectric media of constants ei,e 2 yields 8 



n=l 

This reduces to (106) once one of the slabs is a conductor, i.e Ai = — 1. We have 
therefore provided a derivation of this formula when the conductor is described as 
a statistical system of fluctuating charges in thermal equilibrium. It is interesting 
to note that thermal fluctuations in one of the slabs suffice to generate the correct 
asymptotic value of the force. 
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8 One generalises formulae (5.2), (5.3) and finally (5.5) of [6] starting from (2.4) by keeping ei 
and £2 different. 



23 



A Slab of finite thickness 



The analysis of sections 2-5 applies to infinitely thick slabs. It is interesting to check 
that the asymptotic behaviour of the force does not depend on the slab thickness. 
We consider now that the slab has finite thickness a < oo, while, for simplicity, 
we keep the slab A B infinitely thick. The setting of sections 2 and 3 remains the 
same with the x— integration on limited to the interval —a < x < 0. We then 
follow the same route as in section 4 and 5 by first considering the equation (36) for 
the piecewise-flat profile 

1 f° 

k(x) = 0, x < —a, k(x) = — dx k(x) = ka, —a < x < 

a J -a 

R(x) = 0, < x < d, k(x) = k b , x > d (108) 

It is convenient to choose k(x) equal to the average of k(x) in A^, since we expect 
the latter to be close to its mean value at weak-coupling. 

Solving equation (36) with x' fixed, continuity conditions of x i— > <p(x,x',\t) and 
of x i — ^ d x ip(x,x',k) at x = —a,x = 0,x = d as well with ]im x _>± 00 ip(x,x','k) = 
yields 

{(Paa(x, x', k), — a < x, x' < 

(p AB (x,x' - d, k), -a<x<0<d<x' (109) 

<Pbb{x — d, x' — d, k), d<x,x' 



where 



<PAa(x, X , k) = ^ (kA+k)e k A"a 1 -(kA-k)c- k Ao-a 2 

(110) 

8nk [(k A + k)e kAa e kAX + (k A - k)e- kAa e~ kAX ] e~ kBX ' 
VAB (x,x,k) = — - k)ekAa ^ _ — _ k)e _ kAaa2 (HI) 

WD J T T > l<) - — ( p -k B \x'-x\ , -k B \x'+x\ {k A +k)e k A- aA -( kA -k) c - k A- ai \ , . 

ip B B{x,x,K)- ye +e {kA+k)c k A a ai _ {kA _ k)e - kA a a2 J {lit) 



and 



<3\ = (k A +k)(k B +k)c kd -(k A -k){k B -k)c- kd <J 2 = (k A -k)(k B +k)c kd -(k A +k)(k B -k)c- kd 

(113) 

0"3 = (k A +k)(k B -k)c kd -(k A -k)(k B +k)c- kd (74 = (k A -k)(k B ~k)c kd -(k A +k)(k B +k) e - kd 

k A = \Jk 2 + k a 2 k B = \Jk 2 + k b 2 . 

One deduces also from the differential equations that for any a > both (p and $ 
verify the charge sum rule (34). As a — > oo, formulae reduce to those obtained in 
section 4 for two semi-infinite plasmas. 
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The main observation to be made on this explicit result is that it obeys exactly 
the same factorisation property in terms of the scaled variable q = kd as (44) with 
the the same factors (here (p° A (x, 0, 0) corresponds to the single plasma A with finite 
thickness). The rest of the analysis is the same as in section 4, with the difference 
that the above solution verifies a— dependent bounds in place of (45)- (47), namely 
(appendix B) 

47T 

\<f(x, x', k)| < ^(x, x') < — coth/ta, nd > 1, (114) 

K 

where ^(x, x') is defined piecewise from <p AA , <f AB , etc. as in (109) with 

v. . . 2tt cosh k(cl — \x' — x\) + cosh n(a — \x' + x\) . , 

<p AA (x,x) = — (115) 

AAy ' k smh/ta v ' 

> i i\ 47r cosh«(a - Id) _ / . . 

<Pab(x, x ) = — A r^r — e 116 

K 2 d smh na 

^(x, x ') = *L. [q-^-A + e-^l^'+'l) , K := min{K A , k b } (117) 

The potential $ with structured profiles is related to cp by the integral equation (38) 
and the estimate (49) of the lemma becomes 

\®(x, x',k)\ < 3— — v? > (x, x'), r(a) := r coth/ta (118) 

1 — r(a) 

where r is defined by (48). To have r(a) < 1 one needs r sufficiently small (weak- 
coupling, see appendix B) and na not too small, i.e. the slab width is larger than 
the typical screening length in the plasma 9 . Then the steps leading to (62) are the 
same as in section 4 and the considerations of section 5 apply as well. The reason 
for the asymptotic force being independent of the slab thickness is clearly displayed 
in expressions (62) and (81): it only depends on the screening cloud associated to 
charges located at the inner boundaries of the slabs, and thus is insensitive to charge 
fluctuations elsewhere in the slabs. 



B Bounds for the Debye-Hiickel potential 

In this appendix we present some details of the calculations leading to bounds used 
throughout the paper for the Debye-Hiickel potentials cp and $ and discuss the 
validity of (49) in the weak-coupling regime. 

At first, we show the bound (114), which is a generalisation of (45) to the case 
where plasma A is of finite thickness a. Result (45) is recovered by taking a — > oo. 

9 Notice that the above bounds cannot be uniform in a : ip A (x, x' , 0) diverges as a — > so that 
its integral over [—a, 0] leads to the constant value requested by the charge sum rule. 
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From (113), one has <J\ > cr 2 ; <j\ > /tAKs(e fcd — e kd ), and <t 2 + o"4 < Ci + 03; 
°"3 — °"i < °4 — °"2 < 0. This implies 

1 < (fc A +fc)e fc A"a 3 -(fc A -fc)e- fc A° ( T 4 

— (fc A +fe)c fc ^ a (Ti-(A:A-A:)e- fe A a ( 72 — ' 

which yields the bound (117) for ipsB '■ 

\<Pbb\ < g (e- k ^'^ + e- fc -l^l) < <p% B (x,J). 
To obtain the bound (116) for </?ab, we then note that 

*i - fe| e " 2feAV 2 > (1 " e~ 2feAa ) ai. (119) 

/c 1 1 kd 111 
Thus < < 

' o-, _ >Udl e -2k A a a - 1 _ e -2k A a K K £ Q kd _ Q -kd - 1 _ Q -2k A a ^ 2 ' 

so that 

fe A (a-|a:|) , p -k A (a-\x\) 

< M^.k) < & efcAa I e _ fcAa e" fc - < ^(x.x'). 

Last inequality uses 

<^ 3-^, 0<x<X, k>K>0. (120) 



Finally, (119), the fact that |<r 2 /oi| < 1 and (120) show the bound (115) 
\(p AA (x,x',\s.)\ < 



2 ,j ie - fe A|a;'-a;| _|_ l^jg-feAN'+^l _|_ e -2fcA« Q^jgfcAlx'-xl _|_ a ^ e k A \x'+x\^ 



< ^(i,a/). 
Proof of the lemma 

To proof the bound (49) of $(x, x', k), we proceed as follows. By (38) we can develop 
x',k) as a perturbation series w.r.t. x', k), whose n th term reads 

(471) / dsi " ' ds ™ M ( Sl ) ' ' ' M ( S ™M X > s i> k M s i» s 2, k) • • • (p(s n , x', k). (121) 

This term is bounded by r n tp > (x, x'), where r is given by (48). Indeed, according to 
(45), tp(x,x',\s) is bounded by ip > (x,x'), which itself satisfies 

47T 

ip > (x,s)ip > (s,x') < — v? > (x,x / ), Ws,x,x'. (122) 

K 



26 



Consequently, if r < 1, the series is absolutely convergent and the lemma holds. 
Inequality (122) is proven using (46), (47) and verifying it for each case. As an 
example, 

V > A (x, s) V > AB (s, x') = — 4^e~ K '*'' {e'^ x - s l+M> + e'^^) (123) 

K K (1 

4:71 

< —^> B (x,x'), (124) 

because |x±s| + \s\ > \x\. Some of the majorations leading to (122) assume nd > 1. 
If plasma A is finitely-thick, (122) generalizes to 

47T 

^(x, s)ip > (s, x 1 ) < — coth(/ta) (^ > (x, x'), Vx,x',s (125) 

K 

and we obtain the bound (118) for $. 



Profiles in weakly-coupled plasmas 

The parameter r occurring in the bound (49) can be chosen small enough in the 
weak-coupling regime, defined by T = ^j3e 2 K < 1 (e is a typical charge of the 
system). Indeed, to estimate r, the deviations of the density profiles to their bulk 
values need to be known. In the simplest case of a semi-infinite charge-symmetric 
plasma 10 in the weak-coupling regime, Jancovici [19] finds 

M 7 *)-P, 7 ^ r<<i _ 

PA-y 2 

where x is integrable. Integrating over x < and forming r according to (48) shows 
that r is proportional to T. This also holds if the semi-infinite plasma is not charge- 
symmetric (see [15], sect. 5). For our two-plasmas system, r will be less than 1 
provided T is small and d is large. 



C Decay of Mayer graphs at large slab separation 

We consider prototype graphs constituted of bonds (76) labelled by the indices 
AA, AB, BA, BB according to the respective location of the variables x, x' in slab 
A or in slab B. In view of (29) and after the changes ltd = q, x' — + x 1 — d, the 
contribution to the force of a graph Hab with first root point in and second root 
point in A B is 

-i /»0 /»oo /»oo _ 

(ff AB = ~ T2 dx / dx7 dgge-^l— 'le-^e 7 e y p^( 7 x)p i? (7 , x')7r^(7x,7'x / , 

« J~ooJo JO 7j7 , 

' (126) 

10 A charge-symmetric plasma has two species with opposite charges and same bulk densities. 



27 



A graph Tl AB having L bonds of T AB or T BA type, written in Fourier space with 
respect to the y-variables, is of the general form 

Tl AB ( 1 x, 1 'x>,k) = ^- fdlp{\)--- Jdmp(m) J\[J^T [AB] {£ 3 ) 

(127) 

Here m is the number of internal points, % = (7,, Xj) and J di stands for summation 
over particle species and integration over x^. T(k) stands either for F(k) or F R (k) 
and we have omitted the species and x, x' dependencies from the notation 11 . The 
product of m + 1 ^-functions expresses the conservation of wave numbers at the 
m internal points plus a relation that fixes the sum of ingoing (or outgoing) wave 
numbers to k, as a result of y— translation invariance. These constitute m+1 linear 
equations between wave numbers from the sets {£j}, {k a }, {k b }, which imply C 
independent relations involving only £ variables. Depending on the topology of the 
graph, 1 < C < L. Consider for instance the graph constituted by a single chain of 
bonds F AB (ti)F§ B (k b )F BA (e 2 )F% A (k a )FBA(i3) with L = 3 12 . The conservation laws 
k = £1 = k a = £ 2 = k b = £ 3 imply the independent relations k = £1, £\ = £ 2 , £2 = £3 
between the £ variables, thus C — 3. Consider now the graph constituted of two 
parallel chains Ff 4 (k al )F AB (^ 1 )F| B (k 61 ) and F^ A (k a2 )F AB (£ 2 )F§ B (k b2 ) with L — 2. 
The conservation laws are k = k al + k a2 , k al = £ 1 = k 61 , k a2 = £ 2 = k b2 implying 
the single relation k = £ 1 + £ 2) thus (7 = 1. 

Then we perform the integrations in (127) in the following order. We first carry C 
integrals on the S functions corresponding to the above relations between £ variables: 
this expresses the £j variables j = 1, . . . , L in the integrand as linear combinations 
of the remaining L — C ones, say lc+i, ■ ■ ■ ,£l- We evaluate now Tl AB at k = q/d 
and change the variables £j = Qj/d, j = C + 1, . . . ,L: the Jacobian provides a 
factor d~ 2(yL ~ c \ As d — > 00 the k a , k^ and q integrals factorize. Indeed in T AA or 
J~bb bonds the q dependences occur in the form q/d — > 0, d — > 00 whereas for Tab 
or T B a bonds we use the asymptotic form (60). This yields a factor d~ n ^ AB ^ if the 
number of F[ A b] bonds is n[AB] and d~ in ^ AB ^ if the number of F^ B ^ bonds is n R [A5] 
(for the latter, see at the end of this appendix). Then the k a and k^ integrals refer 
to products of Taa and T BB as in single semi-infinite plasmas and the q integrals 
are carried on product of functions q/ sinhg. For the above examples the q integrals 
are f dq(q/ sinhg) 3 and f dqi f dq 2 (qi/ sinhg 1 )(|q 1 — q 2 \/ sinh |q x — q 2 |). As a final 
result a graph decays as 

d - 2{ L-C) d -nlAB] d -4n*lAB]^ & ^ ^ (12g) 

lx In (127), T[ab] designates either a Tab or a Tba bond. Notice that TBA{i,j, k) = Tab{], i, k). 
12 The F R bonds are needed because of the excluded convolution rule between F bonds. 
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times a factor of order one resulting of the above integrals. It is clear that the 
minimal decay d~ l is obtained when there is only one Fab bond. 

It remains to examine the decay of a F^ B ^ bond, which reads in Fourier form, 



according to (21) 



7 K 
AB 



(k) = /« 



dy e~ lky [exp(-/5 ei e 2 $ AB (y)) - 1 + Pe 1 e 2 ^ AB (y)} . 



(129) 



The x variables are omitted and \x± — x 2 \ is large enough so that the short range 
part of the potential does not contribute. In view of (60), ab (y) has the asymptotic 
form 



$AB(y) 



l 



dq 



d 2 J (2tt) 



exp [q • y) & ab (§) 



!/(*) with /(y) = $«|- 



dq 

(2^ 



47r sinh g 



(130) 



Hence substituting (130) in (129) and expanding for large d after the change of 
variable u = y/d gives 



FAB \d 



~d 2 /due- iq u 



«P I -£/<»)) -l + |/<«) 



1 



due- iqu (/(u)) 2 = C 



(131) 
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